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Pure-jump processes have been increasingly popular in model¬ 
ing high-frequency financial data, partially due to their versatility 
and flexibility. In the meantime, several statistical tests have been 
proposed in the literature to check the validity of using pure-jump 
models. However, these tests suffer from several drawbacks, such as 
requiring rather stringent conditions and having slow rates of conver¬ 
gence. In this paper, we propose a different test to check whether the 
underlying process of high-frequency data can be modeled by a pure- 
jump process. The new test is based on the realized characteristic 
function, and enjoys a much faster convergence rate of order 0(«V 2 ) 

(where n is the sample size) versus the usual o(n 1 ' 4 ) available for ex¬ 
isting tests; it is applicable much more generally than previous tests; 
for example, it is robust to jumps of infinite variation and flexible 
modeling of the diffusion component. Simulation studies justify our 
findings and the test is also applied to some real high-frequency fi¬ 
nancial data. 

1. Introduction. Ito’s semimartingales are widely used in modeling the 
log prices of an asset since they fit many stylized features of asset returns, 
and in option pricing due to absence of arbitrage in efficient market. Mathe¬ 
matically, they consist of two parts: a continuous local martingale term and 
a pure-jump process with both big and small jumps. Ito’s semimartingale 
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with a continuous local martingale is in common use in the literature, for 
example, the Black and Scholes (1973) model (geometric Brownan motion), 
the Merton (1976) model and Kou (2002) model (geometric Brownan motion 
plus finitely many jumps). 

On the other hand, in recent years pure-jump processes have also been 
accepted as an alternative model for log price processes or even the la¬ 
tent spot volatility process to the classic models mentioned earlier; see, for 
example, Todorov and Tauchen (2010, 2014) and references therein. The 
idea behind the pure-jump modeling is that small jumps can eliminate the 
need for a continuous martingale. Pure-jump models are also very flexi¬ 
ble. They include the normal inverse Gaussian [Rydberg (1997); Barndorff- 
Nielsen (1997, 1998)], the variance gamma [Madan, Carr and Chang (1998)], 
the CGMY model of Carr et al. (2003b), the time-changed Levy mdoels of 
Carr et al. (2003a), the non-Gaussian Ornstein-Uhlenbeck-based models of 
Barndorff-Nielsen and Shephard (2001) and the Levy-driven continuous-time 
moving average (CARMA) models of Brockwell (2001) for the stochastic 
volatility. Pure-jump models have been extensively used for general option 
pricing [Huang and Wu (2004); Broadie and Detemple (2004); Levendorskii 
(2004); Schoutens (2006); Ivanov (2007)] and for foreign exchange option 
pricing [Huang and Hung (2005); Daal and Madan (2005); Carr and Wu 
(2007)]. Other applications of pure-jump models include reliability theory 
[Drosen (1986)], insurance valuation [Ballotta (2005)] and hanancial equi¬ 
librium analysis [Madan (2006)]. 

Statistically, this forces us to reconsider the necessity of including the lo¬ 
cal martingale part driven by Brownian motion in modeling high-frequency 
data. This begs the following question: “Is it sufficient to model high fre¬ 
quency data by pure-jump process alone,” or equivalently, “is it necessary 
to add a Brownian force underlying the high frequency data?” The answer 
to this question serves as a model selection purpose. For more motivation 
and explanation, we refer to Ait-Sahalia and Jacod (2010) and Jing, Kong 
and Liu (2012). 

For ease of presentation, let Xt be a semimartingale defined on some 
filtered probability space (DjJ 7 , P), 

X t = X 0 + f\ s ds+ f a s dW s + Xf, 

Jo Jo 

where Xq is the initial value, b s ds is the drift term with b s being the time- 
varying drift coefficient which is an optional and cadlag process, /[[ a s dW s 
is a continuous local martingale with a s being an adapted process and W s 
a standard Brownian motion and the last term is a pure-jump component 
with the jump activity index (3 defined by 

P = inf{r; £ \A a X\ r 

^ 0<s<T 


( 1 . 1 ) 


< oo 
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where A S X = X s — X s _-, see Ai't-Sahalia and Jacod (2009) and Jing et al. 
(2012). Then the above question can be formulated as a hypothesis testing 
problem as 



( 1 . 2 ) 


where T is the time span of the high-frequency data. 

The testing problem (1.2) was studied by several authors. For instance, 
Cont and Mancini (2007), Ait-Shalia and Jacod (2010) used threshold power 
variation to construct their test statistics. However, there are two main 
drawbacks with the threshold power variation method: 

• First, their tests require that X d be of finite variation, which rules out 
many interesting models in finance since empirical evidences in some real 
data analysis show that the jumps are of infinite variation; see, for exam¬ 
ple, Ai't-Sahalia and Jacod (2009) and Zhao and Wu (2009). 

• Second, their tests are not very powerful, even when /3 (0 < (3 < 2) is close 
to 0. This is rather counterintuitive since probabilistically the smaller the 
value of /3 is, the farther X d is from a continuous semimartingale. 

Interestingly, Todorov and Tauchen (2011) invented a test based on the 
point estimator of the JAI cleverly constructed as the smallest power for 
which the realized power variation (without thresholding) does not explode. 
Surprisingly, a test based on this estimator for the presence of Brownian 
motion has the property that it has more power for lower level of activity. 
However, since it is from realized power variation, once more, one has to 
assume that X d is of finite variation when Ct = Jq ds does not vanish 
in order to have available central limit theorem. It is also worth noticing 
that Todorov and Tauchen (2014) test for presence of Brownian motion by 
checking whether “devolatilized” truncated returns are i.i.d. normal assum¬ 
ing finite activity jumps present in the underlying log price processes. 

Testing the existence of a nonvanishing continuous local martingale is 
challenging when the jumps are of infinite variation. Jing, Kong and Liu 
(2012) used the number of small increments to propose a test, which mit¬ 
igates the above-mentioned difficulties, and can handle jumps of infinite 
variation. However, it still has the following deficiencies: 

• First, the local volatility model is too restrictive. For example, it does not 
even cover the Heston model under Hq. 

• Second, the spot volatility of the continuous component is assumed to be 
positive almost everywhere in time t. So if Hq is rejected, it is quite possi¬ 
ble that the continuous component vanishes only in certain subintervals, 
but is still present in other subintervals; see the simulation in Section 4 
for more illustration. 
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In this paper, we develop a novel test to (1.2) to overcome the difficulties 
encountered in previous approaches. The convergence rate of our new test 
under Hq is of order nA 1 / 2 when the jump component is of infinite variation, 
which is faster than that of all existing tests. The idea of the test is based on 
the realized characteristic function, which was introduced in Todorov and 
Tauchen (2012) to investigate the distributional property of volatilities at 
different time points; see also Todorov, Tauchen and Grynkiv (2011) and 
Jacod and Todorov (2014). With observable i.i.d. increments of a class of 
Levy process with either finite activity or infinite activity jumps, Chen, De- 
laigle and Hall (2010) proposed a regression method based on the empirical 
characteristic function to estimate the parameters of the drift, scale, stable 
index and the distribution of the jump size of a compound Poisson process, 
while in our paper, we assume a flexible Ito semimartingale with stochastic 
volatility and stochastic coefficient of jump measures, and assume that the 
time lag of successive observations shrinks to 0 (high-frequency data) rather 
than fixed, as implicitly assumed in Chen, Delaigle and Hall (2010). How¬ 
ever, we remark that direct application of the realized characteristic function 
does not work in testing ( 1 . 2 ), and some other novel statistical techniques 
are needed. 

The paper is organized as follows. In Section 2, we give some assumptions 
and introduce our test statistics. Main results are presented in Section 3. 
Section 4 gives some simulation studies and real data analysis. The main 
proofs are postponed to the Appendix, and the proofs of some lemmas are 
provided in the supplementary material [Kong, Liu and Jing (2015)] to this 
paper. 

Throughout the paper, we assume that the available data set is { X ti ; 0 < 
i <n} which is discretely sampled from X, and is equally spaced in the fixed 
interval [0, T ], that is , U = iA n with A n = T/n for 0 < i < n. Denote the jth. 
one-step increment by 

A]X = X tj -X tj _ 1 , 1 <j<n. 

2. Methodology. The key idea behind our test statistic is that the char¬ 
acteristic function of the increments of the Ito’s semimartingale is dominated 
by the continuous local martingale part. 

For illustration, let us take the following simple example: 

X t = aW t + 1 Y t , 

where a > 0 is a constant spot volatility, 7 is some constant and Yf is a sym¬ 
metric /3-stable process. Then the logarithm of the characteristic function 
is 

( 2 . 1 ) log <p n (u) = log £[ e v /rT “A»X/g^] = _ 


|7|VAi-W a . 
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As A n —>• 0, the last term in ( 2 . 1 ) induced by the jump part decreases at a 
rate of An . Note that when /3 < 1 (i.e., 1 ) is of finite variation), in the 
context of estimating a (or its functionals), the bias caused by the jump 
part is of negligible size o( A n ). This implies that an estimator of at (or its 
functionals) for a general semimartingale based on the characteristic function 
would very likely be robust to jumps of finite variation, which is confirmed 
in Todorov and Tauchen (2012) and Jacod and Todorov (2014). On the 
other hand, the problem becomes more challenging when (3 > 1 since the 
last term in (2.1) is no longer a negligible bias term. In testing (1.2), under 
Hq, the right-hand side of (2.1) is a nonvanishing constant while under Hi 
it is almost zero. This is a major feature we will explore later to differentiate 
the null and the alternative hypotheses. 

We shall now introduce our test statistic. To start with, we split the data 
into m n nonoverlapping blocks with each block length equal to 2v n consisting 
of 2k n intervals of length A n , where k n is some integer depending on n. 
Motivated by (2.1), and in view of X f+S — X t ~ at(Wt+ 8 — IT)) + (l)j ) 5 — 

Y t + ) + 7 t T (Y^s — Y t ~) where Y± are two independent “stable like” Levy 
processes and are two cadlag processes that will be specified later in 
Assumption 3.1. When s is close to 0, we can estimate a'ij Vn (0 < j < m n — 1) 
locally by 

(2.2) c°(u) = Xlo g ( < L SWV-L), 
where 

-j kri 1 

(2.3) L$(«) = —— ^ cos(u(A- kn+2l+1 X - A% jkn+2l X)/A V 2 ). 

n 1 1=1 

Summing over Cj(u) for all j < m n and properly normalizing it, one easily 
gets an estimator of the integrated volatility process, 



Jacod and Todorov (2014) introduced a bias-corrected estimator of Ct as 

[t/(2v n )\-l 

(2.4) Co(u n ) = 2v n Y ( c % u n) - 271 - 7 -(smh(^c°(u n ))) 2 J, 

J= Q V u n\ K n — i-) J 

and further showed that 

(2.5) CoK) = C t + A 0 (u n )t + O p ( A}/ 2 ), 
where 

Ao(u)t = 2u /3_2 A) 1 _/3 / 2 f a s ds 

Jo 
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with a s = x{P){\lt\ P + 1 7 S and x(P) = / 0 °° V P sin ydy. Then a natural 
test statistic which can differentiate the null and alternative hypotheses is 

, = Cp(2u n ) - Co(u n ) ( 0, on {C T > 0}, 

n_ C 0 (u n ) P l2 M -K0, on {Ct = 0}. 

The problem with T.'' n is that no central limit theorem is available as (3 > 1, so 
that one cannot find the rejection region when jumps are of infinite variation. 
We will fix this problem with some manipulations to T' n below. 

To do this, we replace Cq(u) by a similarly defined quantity. Let the cj’s 

and Ci(it) be similarly defined as the Cj’s and Co(u) with A 2 jk n + 2 i+\X ~ 
A 2 jk n + 2 i X replaced by A% jkn+2l X - A% jhn+7l -\ x , tor 1 = 1 ,... X -l.A 
seemingly better test statistic is then 

tr, r \ (C 0 (2u n ) — Ci(u n )) — (Co(2u n ) — Co(u n )) C 0 (u n ) - Ci{u n ) 

( 2 . 6 ) =---=---, 

CiK) C^Un) 

which works under Hq because the numerator is equal to 

[(0,(210 -Ct- Ao(2u n )?) - (CiK) -Ct- A 0 (u n )?)\ 

(2.7) -[(C 0 (2u n ) - C t - A 0 (2u n )?) - (C 0 (u n ) -C t - A 0 (u n )?)\ 


= O p (A 1 J 2 )-o p (A 1 n / 2 ). 


1 /2 ^ ^ 

The second term in (2.7) is o p (A n ' ) since Co(2u) and Co (it) are calculated in 
the same way, except for using different arguments, and are asymptotically 
perfectly correlated as it = u n —>• 0; see also (a) in Theorem 1 of Jacod and 
Todorov (2014). However, the first term in (2.7) is O p (An ) since Ci(it n ) 
uses the data points one grid after those in Co(2it n ), which decreases the 
overlap of the data and hence has lower dependency between the terms with 
argument 2u n and u n ; see Theorem 3.2 below. 

Although T* / An is tight under H$ , it can be close to zero with a large 
probability under H\ since the signal in the numerator is swept away in the 
bias correction. This causes difficulty in successfully detecting pure-jump 
processes under H\ and hence results in a low power. This difficulty can be 
remedied by adding a bias of order o( A n ) onto the numerator of T*. 

Our final test statistic is 


( 2 . 8 ) 


T n 


C 0 (u n ) - Ci(u n ) - y n A 


1/2 

n 


Cl (Un) 


where is some chosen constant satisfying —>• 0 of which the explicit 

form will be given in Section 3.3. It can be shown that 


rp / A 1/2 / = Opi 1 ), on {Ct > 0 }, 

n/ n X P -oo on {Ct = 0 }. 


(2.9) 
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This means that T n / A 


1/2 

n 


can be used to differentiate Hq and H\. 


3. Main results. 


3.1. Model assumptions. We need the following assumptions. 


Assumption 3.1. 

x t = [ 7 f-dY s + + [ 7 ~_dY~+ [ f <5(s, z)p(ds, dz), 

Jo Jo Jo Jr 

where Y + and Y~ are two independent Levy processes with positive jumps 
and Levy triplet equal to (0,0, F±), 7 ± are two cadlag adapted processes and 
p is a Poisson random measure on R + x R with intensity q(dt, dx) = dt®dx. 
We assume further that, for some (3 > 1 > r, the Levy measure satisfies 



^((x.oo)) 


1 


<9{x), 


xG ( 0 , 1 ], 


with g{x) a decreasing function s.t. x r 1 g(x) dx < 00 , and 1 5(t, x)\ r A 1 < 
J{x) with J{x) Lebesgue integrable on R. 


Assumption 3.2. at is an Ito semimartingale of the form 

ft rt ft 

a t = a 0 + 


Jo 


bg ds + 


H°dW s + 


H' s a dW' s 


to 


/ 0 


+ f f S a (s,x)(p - q)(ds,dx) 

JO ^{|£ £T (s,rr)|<l} 


+ 


ft 

Jo J{ I 


{| <S cr (s,a;) | >1} 


^(s, x)p(ds, dx), 


where all the integrands are optional processes satisfying the integrable con¬ 
dition in Ito’s sense, and q is the compensator of p. Assume that W and 
W' are two independent Brownian motions that are further independent of 
(p,Y+,Y~). 


Assumption 3.3. We have a sequence r n of stopping times increasing 
to infinity, a sequence a n of numbers and a nonnegative Lebsgue-integrable 
function J on R, such that the processes b, H a , ^ are cadlag adapted, 
the coefficients 5, 5 a are predictable, the processes b a , H ,cr are progressively 
measurable and 

t<r n => \S(t,z)\' A 1 < a n J(z), |<5 CT (t,z )| 2 A 1 < a n J(z), 
t < T n , V — b, b°,H°, H' a , 7 ± =>\V t \< a n , 

V = b,H *, 7 * =7 | E(V {t+s)ATn - V tATn \Ft)\+E(\V (t+ s)A Tn ~ VWJ Vt) < 
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Assumption 3.1 is the same as the Assumption (A) given in Jacod and 
Todorov (2014). It essentially states that X d can be decomposed into two 
components: active and less active jumps. Here, the first two components are 
the stable-like jumps assumed to have the jump activity index /3 > 1. (This 
can be extended to cover the case for r < f3 < 1 with extra efforts and possibly 
more stringent conditions. However, if we have a priori that /3 < 1, more 
straightforward tests will be possible.) Another reason we restrict attention 
to /3 > 1 is because this is more interesting and challenging statistically. The 
last term consists of jumps with finite variation (but possibly of infinite 
activity) which is expected to disappear in a limiting sense as inspired by 
the finding following (2.1). In A'it-Sahalia and Jacod (2010), it is assumed 
that (3 < 1 since otherwise no asymptotic distribution theory could be used 
under H$ to calculate the rejection region. 

Assumption 3.2 is a standard assumption in the literature which allows 
for the “leverage” effect due to the common driving forces in X and a. In 
Assumption 3.2, the jumps of at are assumed, without restriction, to be 
driven by the same Poisson measure as X. 

Assumption 3.3 is the same as the Assumption (B) in Jacod and Todorov 
(2014) and a rather general assumption which is satisfied by the multifactor 
stochastic volatility models that are widely used in financial econometrics, 
for example, the popular affine jump diffusion models in Duffie, Pan and 
Singleton (2000). Assumptions 3.2 and 3.3 admit a rather general Ito semi¬ 
martingale as the continuous part under Hq. As a comparison, Jing, Kong 
and Liu (2012) require that the volatility be of form a(Xt), a smooth func¬ 
tion of Xt bounded away from 0. Hence our assumptions on the continuous 
component is far less restrictive than that in Jing, Kong and Liu (2012). 

3.2. Main theorems. We first state a central limit theorem for the joint 
distribution of (Co(u n ),Ci(u n )). 


k n X 


Theorem 3.1. Suppose k n , u n , 7 n and A n satisfy 

k n X l J 2 ~ £ -A OO, U n -A 0, sup 

n 'LL, 

(3.1) 

7 n —•► 0 , 

for any e > 0. Let c s = a 2 . Then on the set {Ct > 0} we have 

1 ( C 0 (u n ) - A 0 (u n )? - C t ' 

A i/2 {c l{Un )-A 0 (u n )^-C tj 

f C s dW 8 

Jo 


1/2 


< OO, 


( 


\ 


K j\ s d(^W s + V3/2wA ) 


(3.2) 
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where W and W^ are two mutually independent standard Brownian motions 
defined on an extension of the original probability space and are further 
independent of T, and C s stands for stable convergence. 

In Theorem 1 of Jacod and Todorov (2014), a similar multivariate central 
limit theorem related to the bias corrected estimator of Ct in (2.4) with 
distinct arguments was obtained. While in (3.10) and (3.11) of Theorem 1 
of their paper, the vector of component estimators with distinct multiples 
of u n are formed by using the same way of aggregating the high-frequency 
data, Theorem 3.1 in our paper considers a bivariate central limit theorem 
for (Co(u n ), C\(u n )), with Co(u n ) collecting the high-frequency data one 
lag after Ci(u n ). By simple application of Theorem 3.1 and the continuous 
mapping theorem, we soon have the following null distribution of T n . 


Theorem 3.2. Under the conditions in Theorem 3.1, we have in re¬ 
striction to {Ct > 0}, 


a~ 1/2 t„ 


G, 


where Gt is a centered Gaussian process with conditional variance Kt = 



It follows from Theorem 3.2 that the convergence rate of T n is of order 
A, 1 / 2 , in contrast to A^ 4 in Jing, Kong and Liu (2012), where w > 
(3 — 1/2 is some constant (practically w is taken as 3/2 since (3 is usually 

of /q 

unknown) or Vn ' in Ait-Sahalia and Jacod (2010), where (3' < 1 and v n 
satisfies 


v n /AO —)• 0, v n /A(/ + —»■ oo, 0 < p- < p+ < 1/2. 


Theorem 3.2 is not directly applicable in determining the rejection region 
since the conditional variance is unknown. The denominator of the condi¬ 
tional variance can be consistently estimated by (Ci(u n )) 2 , thanks to (2.5). 
Inspired by the construction of Ck(u ) (k = 0,1), we use the following linear 
combination of sample variances to estimate the integral in the numerator 
of kt- Define 

(3-3) i n = + In,l)i 

where 


(3.4) 



(sinh(u 2 c^(u n )))\ 2 

uUK-i) )' 


A; = 0,1. 


Now we have the following studentized central limit theorem. 
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Theorem 3.3. Let kt = 4/ n /(Ci(u n )) 2 . Then we have under the con¬ 
ditions in Theorem 3.1, in restriction to {Ct > 0}, 


(3.5) 


T n = 


1 


A 1 / 2 k V 2 
L-\n tv'T' 


Co(u n ) - Ci(u n ) - 7 „A 
27 1/2 A 1/2 

4±n L-±n 


1/2 


A7(0,l), 


where jV(0,1) is a standard normal random variable independent of F. 


From Theorem 3.3, we can reject Hq ifT n < ~z a where P(AT(0, 1) > z a ) = 
a for a € (0,1). Now we state a result on the convergence rate of T n under 
Hi. 


Theorem 3.4. Suppose Assumptions 3.1 and 3.3 hold, k n A.)/ 2 — > 0, 
k n A l n /2 £ —> oo for any e > 0, sup n k n A l J 2 /uif < oo and u n is bounded. Then 
on the set {Ct = 0, f^a s ds / 0}, we have 


(3.6) c„K) - OK) = o p (u- n 2 M~ mi,+1 - rn + »W 2 - 2 Ai-W) 

and 

(3.7) /„=4„f-‘A 2 < 1 -« 2 > f alds + o r (uf-‘C$'-t>IV). 

Jo 


The following result concerning the size and power performance of the 
test is a straightforward consequence of Theorems 3.3 and 3.4. 

Corollary 3.1. (1) Under the conditions in Theorem 3.1, we have 

P(T n < —z a \{CT / 0}) -A a; 

(2) under the conditions in Theorem 3.f, if 

,,(«»+>-« + u 2- W 2 A «4-l/2 ) _ ^ 

we have P{T n < —z a \CT = 0, a s ds / 0) —>• 1. 

Remark 3.1. Corollary 3.1 shows that our new test achieves asymptotic 
nominal level a and the asymptotic power 1. It follows from the proof of 
Corollary 3.1 that T n goes to —oo with rate O p { 7n(^)^ 2_ ^ //2 ) under H\ 
and conditions in 2. Thus the test becomes more powerful as f3 gets closer to 
0, which will be further confirmed by our simulation studies. This overcomes 
the drawbacks of the test by Ait-Sahalia and Jacod (2010). 
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3.3. Choice of tuning parameters. We now study how to choose tuning 
parameters k n , u n and j n . The major role of k n is to balance the bias and 
variance of Co(u n ) — Ct and Cj (u n ) — Ct- The larger the k n , the smaller the 
bias and the larger the variance. Hence we could choose k n = —c' A n '“ log A n 
for some constant d > 0. 

Now we turn to u n . The rationale for letting u n —> 0 under Hq is to guaran¬ 
tee the convergence in probability in (A. 13). As in Jacod and Todorov (2014), 
we choose u n so that u 2 f 0 c s ds —>• 0 by setting u n = c(log (1/A n )) -1 / 30 x 
where BVt = (7t/2) Y^i=i |A”X|| A” +1 X| * s the bipower varia¬ 
tion, which is a consistent estimator of c s ds. Another advantage of such 
choice of u n is that it would be enlarged under H\, which in turn increases 
the power as is seen from Corollary 3.1 and Remark 3.1. Choosing an optimal 
c is quite hard. In order not to incur much approximation error in (A. 13), we 
suggest to choose small c when n is moderate, say c = 0.18. Simulation stud¬ 
ies where the data is generated from a fitted model (no guarantee of good 
fitting accuracy) assuming Hq given in Jacod and Todorov (2014) show that 
choosing c around 0.18 would work well. 

Finally, we look at y n . On the one hand, should be close to 0 under 
Hq in order not to produce a big bias for T n ; on the other hand, 7 n should 
converge to 0 with a rate of u~ 2 Af/ 2- ^ 2 ^ 4-1-7 ^ + Un" 2 2 A^ 2 ”^ 4 so that the 
test has good power. This is easily achieved by setting 7 n = c*/log (u^/ A n ) 
when u n is determined by the aforementioned method. To be conservative, 
one can choose small c* when n is moderate, say c* = 0.2. 

4. Numerical experiments. 

4.1. Simulation studies. In this section, we conduct simulation studies 
to check the performance of the new test and make comparisons with the 
test given in Jing, Kong and Liu (2012). We first consider the performance 
on control of type I error probability. As in Jacod and Todorov (2014), 
we generate simulation data for 5000 times from the following stochastic 
volatility model: 

(4.1) Xt = X o + f y/~cfdW s + 0.51t, 0<t<T, 

Jo 

(4.2) c t = CQ+ [ 0.03(1.0 — c s ) ds + 0.15 [ y/T s dW' s , 

Jo Jo 

for 0 < t < 3T/4 and ct = 0 if 3T/4 <t<T. In order to incorporate the 
leverage effect, we set corr (dW,dW') = —0.5. The parameters in the volatil¬ 
ity dynamic are specified by fitting actual financial data in the same refer¬ 
ence paper. The volatility ct is a square root diffusion process which is widely 
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Table 1 

Empirical sizes and the empirical powers of the new test; the nominal level is 5%; 
(n = 1170, k n = 50); (n = 2340, k n = 78); (n = 4680, k n = 100) 


p 


Empirical sizes 


Empirical power 

n = 1170 

n = 2340 

n = 4680 

n = 1170 

n = 2340 

n = 4680 

1.0 

0.0610 

0.0586 

0.0574 

0.9988 

0.9998 

1.0000 

1.1 

0.0616 

0.0624 

0.0610 

0.9984 

0.9990 

1.0000 

1.2 

0.0640 

0.0635 

0.0634 

0.9936 

0.9986 

0.9996 

1.3 

0.0604 

0.0601 

0.0608 

0.9596 

0.9948 

0.9986 

1.4 

0.0522 

0.0616 

0.0616 

0.6508 

0.8414 

0.9650 

1.5 

0.0566 

0.0624 

0.0610 

0.2902 

0.3810 

0.5290 

1.6 

0.0612 

0.0514 

0.0524 

0.1328 

0.1698 

0.2138 

1.7 

0.0594 

0.0624 

0.0554 

0.0942 

0.1068 

0.1208 

1.8 

0.0578 

0.0550 

0.0594 

0.0776 

0.0804 

0.0804 

1.9 

0.0572 

0.0568 

0.0558 

0.0748 

0.0790 

0.0728 


used in financial applications. We tuned k n , u n and as in Section 3.3 with 
c = 0.18 and c* = 0.2. We consider n = 1170,2340,4680 which corresponds 
to sample the data per 20,10,5 seconds, respectively. In the simulation, we 
let T be one day consisting of 6.5 trading hours. 

Table 1 displays the empirical sizes of the new test. Clearly, they are 
slightly higher than the nominal level but acceptable across the board due 
to the small bias added artificially. Figure 1 gives the QQ-plot of the test 
statistics for n = 2340 and /3 = 1.2,1.5, showing that the normal approxima¬ 
tion works well. 

For comparison, we choose f3 = 1.2 and n = 2340 and carry out the test 
given in Jing, Kong and Liu (2012), referred to as JKL’s test below. No 
comparisons will be made with the test given in Ait-Sahalia and Jacod 
(2010) (AJ’s test), since it was outperformed by the JKL’s test in extensive 
simulation studies given in Jing, Kong and Liu (2012). Table 2 lists the 
empirical sizes of JKL’s test where 5* is a tuning parameter determining 
how many small increments are used to compute the test statistics. Clearly, 
the JKL’s test is too liberal since the type I error probabilities are out 
of control, showing that the JKL’s test fails when the continuous process 
vanishes in some subintervals. The reason for the failure is that the JKL’s 
test statistic has a nonnegligible bias, even for large enough n. 

It seems that choosing S* small would have satisfactory control of type I 
error. However, when 6* is small, the normal approximation is actually no 
longer reliable. For 5* = 0.05, there are roughly 5 small increments (effective 
data) used in calculation of the test statistics, which affects the accuracy of 
the normal approximation. Figure 2 gives the QQ-plot for the test statis¬ 
tics given in Jing, Kong and Liu (2012) for d* = 0.05 (left panel), 0.5 (right 







Quantiles of Input Sample Quanties of Input Sample Quantiles of Input Sample 
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QQ Plot of Sample Data versus Standard Normal 



QQ Plot of Sample Data versus Standard Normal 



QQ Plot of Sample Data versus Standard Normal QQ Plot of Sample Data versus Standard Normal 




QQ Plot of Sample Data versus Standard Normal 



QQ Plot of Sample Data versus Standard Normal 



Fig. 1. QQ-plot of the test statistics under Ho for (3 = 1.2 (left panels), 1.5 (right panels); 
from top to bottom, c = 0.15,0.18,0.2; n = 2340. 


panel) when (3 = 1.2 and n = 2340. From the left panel, we see a clear con¬ 
cavity pattern, which implies that the distribution of the test statistic is 
left-skewed, yet the empirical size is 0.07. Apparent improvement in skew- 
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Table 2 

Empirical sizes of JKL’s test; (5 = 1.2, n = 2340; the nominal level is 5%; Empirical 
sizes * stand for the empirical sizes when c s follows the same square root process for 

3/4 T <t <T 


5 * 

0.50 

0.75 

1.00 

1.25 

1.50 

1.75 

Empirical sizes 

0.3032 

0.4442 

0.5916 

0.7312 

0.8532 

0.9402 

Empirical sizes* 

0.0298 

0.0400 

0.0358 

0.0406 

0.0402 

0.0436 


ness could be seen in the right panel for 5* = 0.5 since more effective data 
(roughly 40) were added in calculation of the test statistics. However, we 
see a clear bias in the QQ-plot. 

Next we investigate the power of the new test. We generate the data for 
5000 times from the above model, except that c s = 0. The empirical powers 
for various /3 values are given in Table 1. We make the following observations: 

• the power of the new test decreases as /3 increases since, as /? increases to 
2, the pure-jump process fluctuates more like a Brownian motion; 

• as the sample size increases, the empirical power increases overall, as can 
be expected. 

We also did a sensitivity study to k n when it is chosen in the proposed 
range. In the sensitivity study we take c* = 0.2 and k n = 50,78, c = 0.15 or 
0.2 when n = 2340. The results on both the size and power performance are 
reported in Table 3, where we can see that the empirical sizes and power 
do not change much. We also conducted other sensitivity studies for c ~ 
0.18 and n = 1170 with k n in the corresponding range and reached similar 
conclusions (hence not presented here). 


QQ Plot of Sample Data versus Standard Normal QQ Plot of Sample Data versus Standard Normal 



Fig. 2. QQ-plot of the test statistics given in Jing et al. (2012) for 5* = 0.05 (left panel), 
0.5 (right panel) under Hq when /3 = 1.2; n = 2340. 
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Table 3 

Empirical sizes and the empirical powers of the new test for different pairs of (c, k n ); the 

nominal level is 5%; n=23f0 


p 


Empirical sizes 


Empirical power 


(0.15,50) 

(0.15,78) 

(0.2,78) 

(0.15,50) 

(0.15,78) 

(0.2,78) 

1.0 

0.0630 

0.0604 

0.0634 

0.9986 

0.9994 

0.9992 

1.1 

0.0604 

0.0604 

0.0608 

0.9984 

0.9986 

0.9992 

1.2 

0.0634 

0.0618 

0.0624 

0.9970 

0.9982 

0.9982 

1.3 

0.0558 

0.0592 

0.0638 

0.9842 

0.9896 

0.9968 

1.4 

0.0580 

0.0562 

0.0618 

0.7432 

0.7708 

0.8786 

1.5 

0.0584 

0.0560 

0.0614 

0.3148 

0.3242 

0.4146 

1.6 

0.0576 

0.0613 

0.0608 

0.1670 

0.1498 

0.1814 

1.7 

0.0558 

0.0496 

0.0568 

0.0908 

0.0906 

0.1102 

1.8 

0.0558 

0.0540 

0.0582 

0.0780 

0.0778 

0.0788 

1.9 

0.0542 

0.0544 

0.0566 

0.0702 

0.0702 

0.0744 


4.2. Real data analysis. In this section, we implement our test on some 
real data sets. We first investigate the stock price records of Microsoft 
(MFST) in two trading days, December 1, and 12, 2000, which were also 
included in Jing, Kong and Liu (2012). All data sets are from the TAQ 
database. As in Jing, Kong and Liu (2012), to weaken the possible effect 
from microstructure noise, we sample observations every 1/3 minutes. Fi¬ 
nally, we use logarithms of the sampled prices to calculate the test statistics. 

We set T = 1 (day) consisting of 6.5 hours of trading time. As in the 
simulation studies, we set k n = 50 and 7 n = 0.2/log (u^/ A n ). To be on the 
safe side, let u n take values in the grid points in (0,1] with step length 
equal to 0.01. Figure 3 plots the test statistics against u n for two data sets. 
We see from the figure that for all configurations of u n , the test statistics 
are far lower than —1.645, hence providing significant evidence against the 
existence of a Brownian force. This confirms the empirical results in Jing, 
Kong and Liu (2012) and in the meantime rules out the possibility that 
Brownian force exists in some subintervals. 

Next we implement our test the S&P 500 index data which are sampled 
every 5 minutes during January 4-29, 2010. The tuning parameters are 
used as given above for those two stock data. The observed test statistics 
are plotted against u in the lower panel of Figure 3. We obtain the same 
conclusion that during the specified time period, the underlying log price 
should be modeled by a pure-jump process. 

5. Conclusion and discussion. In this paper, we have developed a new 
test based on the realized characteristic function to check whether the un- 
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U 



Fig. 3. Observed test statistics for the trading date, December 1 (middle panel) and 
December 12 (upper panel), in 2000, and the 5-mins S&P 500 index (lower panel) data 
during January )-29, 2010. The horizontal line has level —1.645. 


derlying process of a high frequency data set can be modeled as a pure-jump 
process, and shown its advantages over existing tests. Here are some future 
problems worth pursing in future research work: 

• The effect of the microstructure noise, in the testing problem (1.2) or even 
in estimating the functionals of the volatility, is unclear and worthy of 
investigation in both theory and practice. Here we could explore the two- 
time-scale technique or multi-time-scale technique [Ait-Sahalia, Mykland 
and Zhang (2005), Zhang (2006)] or the pre-averaging approach [Jacod 
et al. (2009)]. 

• In the present paper, our inference is with the price process. It is of interest 
to make inference on the volatility process which, as recommended in 
Todorov and Tauchen (2014), could be modeled by a pure-jump process. 
The challenge of this problem is that the volatility process is unobservable. 
Studies on this topic is still undergoing. 
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APPENDIX: PROOFS OF MAIN THEOREMS 

This appendix contains the proofs of main theorems. The proofs of Lem¬ 
mas A.4-A.6 as well as some interesting supplemental lemmas are given in 
Kong, Liu and Jing (2015), a supplementary material [Kong, Liu and Jing 
(2015)] to this paper that is not for purpose of publication. By the stan¬ 
dard localization procedure, it is enough to prove the main results under 
the following strengthened assumption. 

Assumption A.l. b, a , 7 + , y - , b a , H a and H' a are bounded. 

Before we prove the theorems, we introduce some notation and give an 
outline of our proof. Let Ut(u) = exp (—u 2 q — 2 An ^ 2 u^a t ) where at = 
x(/3)(| iff + I lt\ P ) with xih = fo°° s[n (y) dy- For ease of notation, 
Uj(u) = U 2 j Vn {u) and sometimes we write Ejr t V s = E(V s \Et) for a stochas¬ 
tic process V*. Let £,k,j( u ) = E k -{u)/U){u) — 1, k = 0,1. Let Q(k,n,t) = {cu, 
maxfcj |^,j(u,w)| < 1/2}. By Lemma 7 of Jacod and Todorov (2014), 

(A.l) P(H c (fc,n,f))^0, 

irrespective of whether the continuous component exists or not. 

A.l. Proof of results under Hq. Assuming the continuous local martin¬ 
gale exists, our proof depends heavily on the following decomposition: 

Cj{u) = c 2jVn + 2uP- 2 A 1 n -P /2 a 2j v 7l - + -^k,j( u ) + r k,j(u), 

(A.2) 

where r^jiu) represents the remaining term which will be shown to be neg¬ 
ligible. By summing up the terms in (A.2) over j, one soon has 

[t/( 2v n )\-l [t/(2v n )]-l 

C k {u ) = ^2 2v n C 2 jv n + ^2 2v n(^~ 2 A 1 ~ P/2 a 2 jv n ) 

j =0 j =0 

[t/( 2u„)]-l 

(A.3) - ^2 2 Vnik,j{u)/u 2 

j =0 

[*/(2«n)]-l ,£2 J u \ 1 x 

+ 2v " ~ _ 1)u 2 ( s mh(u 2 <j( M ))) 2 j + R k (u). 

We will first show that the first and second term converge to some limits, 

1/2 

and the fourth and last term in (A.3) are o p (A n ), while the third term is 

O p (A n ) and converges to a conditionally centered Gaussian random vari¬ 
able stably. This proves the univariate central limit theorem in Theorem 3.1. 
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After that we proceed with the proof of the bivariate central limit theorem 
by investigation into the covariation of those two marginal sequences, which 
ends up with Theorem 3.1. Theorem 3.2 is a consequence of Theorem 3.1 
and the continuous mapping theorem. Theorem 3.3 can be proved by show¬ 
ing that kt is consistent to kt■ In the sequel, K will be a constant that has 
different values at different appearances. 

We now cite three lemmas from Jacod and Todorov (2014), whose proof 
can be found in the same reference paper. Lemma A.l is concerned with 
the first and second term in (A.3), that is, the discretization error terms. 
Lemma A.2 gives the stochastic order of £fcj(a n ), k = 0,1, while Lemma A.3 
shows that the fourth term and the remainder term in (A.3) are asymptot¬ 
ically negligible. 


Lemma A.l [Lemma 8 in Jacod and Todorov (2014)]. Under Assump¬ 
tions 3.1-A.l and assuming (3.1), we have 


(A.4) 


[t/(2v n )\-l 

E 

i=0 


2v n )\ 1 /»£ 

E 2v nC2jVn~ C s ds = O p (u 2 n A)/ 2 ), 

Jo 


[t/( 2v n )]-l 

(A.5) E 2v n{2uP- 2 Al-P /2 a2jV n)-Ao(u n )? = °p(ulAl/ 2 ). 

1=0 

Lemma A.2 [Lemma 14 in Jacod and Todorov (2014)]. Under Assump¬ 
tions 3.1-A.l and assuming (3.1), we have, for k = 0,1, 


(A.6) 

(A.7) 


\E T 2 iv J kd (u n )\<KutAl/ 2 ^ 


Uj(2u n ) + Uj( 0 ) - 2 U 2 (u n ) 


E F 2jv J k ,j(u n ) 2{k n — l)U 2 (u n ) 

and for q> 2, 

(A.8) E T2jVn |&jK)| 9 < K(u 2 n q /k q J 2 + u 4 n v n ) 

where 4> n is some sequence of numbers converging to 0. 


< KuflA)! 2 (\> n , 


Lemma A.3 [Lemma 9 in Jacod and Todorov (2014)]. Under Assump¬ 
tions 3.1-A.l and assuming (3.1), we have R k {u) = Opiu^A^/ 2 ) k = 0,1 and 


(A.9) 


[t/{ 2v„)]-l 

E 2vn 

1=0 


2 < 


1 


{k n - 1 )ul 


(sinh (u 2 n c k Au n ))) 2 ) = o v (u 2 n A l J 2 ). 
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The following lemma provides a formula for the limit of the conditional 
real part of the characteristic function of a linear combination of three suc¬ 
cessive increments. The proof can be found in the supplementary material 
[Kong, Liu and Jing (2015)] to this paper. 


Lemma A.4. Letu* n = |a nj o| V|a n) i| V|a nj 2 |, under Assumptions 3.1-A. 1, 
and assume (3.1) with u* replacing u n , so we have 


^(»-i)A„ cos ( y^ a w 


j=o 


aX 2 / 


1 


ex P i )A , 




1=0 


(A.10) + An 0/2 X(P)j2(\ a n,n^. 1)A J 0 + \ a n,a {i _ 1)An f) 


1=0 

2 


x cos Ai p/2 x'(I3) ^({a n ,i7(t_i) aJ^ + 


1=0 


< KvtfAWfa, 

where {,x = sign {x)\x\^ and x'(/3) = / n °° 1 ~ cos ^'l 


Jo yP 


dy. 


Proof of Theorem 3.1. By Lemmas A.l, A.3 and (A.3), it suffices 
to prove that 


1 /[t/ (2u n )J —1 [t/ (2v n )\—l \ 

A 1 / 2 £ 2'U7 l ^0j(^n)/^rD ^ ^ j(v,n)/u n | 

An \ j =o j=0 / 

converges to the right-hand side of (3.2) stably. By Lemma A.2, we have 

[t/( 2v n )\-l 

^ ] 2 v nE(£k,j( u n)/ u n \J'2jv n ) = °p(. u n ^n )> /c = 0, 1. 

j =0 

Hence it is enough to prove the bivariate central limit theorem with stable 
convergence for the following centered discrete bivariate martingale with 
respect to (J r 2j Vn fj=Q Vn ^~ 1 '- 


2 tv 


A. 


1/2 


(Co ,j( u n)/u n — E{(,0,j (u n )/u n \J-2 j Vn))i 

v j=0 


(A.ll) 
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[t/( 2v„)]—1 

El (Cl ,j( u n)/ u n 
3=0 


E (Cl ,3 ( u n )/ u n I -^2 j v„ ) ) 


Let 


n,0 2Un 

X 3 ~ a 1/2 

L -^n 

n, 1 2"U n 


(Co j ( u n) / u n (Co ,j ( u n ) /\ | J~2jv n ) ) j 


X/ = -i7a(Clj(«n)/«n - ^(Clj(«n)/«„|-? 7 2jt;„)). 

A n 

By Theorem 7.28 in Chapter IX of Jacod and Shiyayev (2003), we only need 
to prove that 


(A.12) < 


sup 

t 


[t/(2u n )]-l 

E 

5=0 

[*/( 2 w «)]—1 - t 

E ^ 2 . r „ n (x”’ fc ) 2 4 J C 2 s ds; 


k = 0 , 1 , 

A: = 0,1, 


3=0 
[t/{ 2v n )]-l 


— /*£ 

E (x?°x- a ) 2 I C 2 S ds ; 

n ^0 


^ 0; 


A: = 0,1, 


3=0 

[t/( 2v n )]-l 

3=0 
[t/(2v n )] — l 

E Ejr 2 j Vn (Xj’ k ( M 2(j+i)v„ ~ m 2jv„)) ~> P 0 ; k = 0 , 1 , 

3=0 


for any square-integrable martingale M. The Hrst equation holds automat¬ 
ically since (Xj’ k )f=o V "^ 1 form a sequence of ^ 2 (j+i)t;„-martingale differ¬ 
ences. 

Now we calculate the conditional variances of the marginal sequences. By 
(3.1), Lemma A.2 and the fact that \Uj(u n ) — e < KAh ^ 2 Un , 

we have 

2 [t/(2^n)] — 1 

EE E ^((CfcbK) - ^(CfejK)|^J) 2 |^J 

anU " J=0 

2 [^/( 2^71 )] 1 

= xAr X - {£[&>„)l^„„]) 2 ) 

J=0 
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4v; 


2 W( 2«n)]-l 




-^{£,k,j iV / n‘')\J~2jv n ) + °p( 1) 


(A.13) = 


l=o 

4v 2 [V(2^)]-i jj.(2u n ) + 1 - 2 Uf(u n ) 

L -- tt 2 i ^ -+ °p( 1 ) 

1=0 


2{k n - 1)A n < j^ n Uf{u n ) 


[t/{ 2v n )]-l 


(1 + °i >(1)) (Uj(2u n ) + 1 — 2U? (u n ))2v n + o p (l) 


U, 


fo( e ~ 4u " Cs + 1 - 2 e~' Zu ^)ds 


3=0 

n [t/{2v n )]-l 

{ + Oj±_» ( e ~4< C 2jvn + 1 _ 2e - 2U « C2 ^n )2u n + 0p( 1) 


1=0 


-,-2 n? < 




+ O, 


Jo 


c 2 s ds , 


where in obtaining the convergence in probability, we used the Taylor ex¬ 
pansion of e 31 when x is near 0. This proves the second equation in (A. 12). 

Next, we are going to check the third equation in (A.12). By Lemma A.2, 
we have 


(A. 14) 


4u 


2 [t/(2v„)]-l 


A n ui 


Ejr 2jv n (£0,jOn) - Ejr 2jvn ^j(u n )) 


1=0 


X (£l ,j( u n) Ej^ 2jVn ^l ,j(u n )) 

4 y 2 [*/ ( 2 ^n )] — 1 

(Ej r 2 j Vn €o,j( u n)€l,j( u r 


A n ui 


1=0 


E J~ 2 jv n €o,l { u n)Ej^ 2jVri £l j {u n )) 


4v. 


2 [* / (2i)n)] — 1 


A n U n 


E^jvn^Oj { u n)£l,j ( u n) + °p(l)- 


1=0 


Now we investigate the summand in (A. 14). Let 
AT 


/ A n . 

> / • 7\ / ‘2jk n -\-2l — fc+1 

CfcU> 0 = COS I U n - 


* - A 21 fen+ 2I- fe * 


1/2 


A. 


Ej r (2jk n +2l-k-l)A n C0S I U n 


/\n y_ An y 

AA 2j/c n +2/-/c+l yv LA 2j/c n +2Z-A;^ Y 


A. 


1/2 
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and 


( A n v A n y\ 

^ 2 jkn+ 2l-k+l ^ 2jk n+ 2l-k \ _ 

k = 0,1. By (6.22) and (6.29) in Jacod and Todorov (2014), we have 

<Ku 4 n A^ 2 <^ + \U 2jv n +(2l—k— 1)A„ (^n) Uj(u n ) |- 
which, together with Lemma A.2 and the property of Ut(u n ), shows that 

/n /(T„ - 

E^jvn £,l,j( u n) | 


(A.15) 


/ (!/(&„-i))E£i 

V Uj{u n ) ) 


s \l E Ki»„K,i( u n)\l E *. 


(i/(fc»-i))Efer 1 ic;o,o-coO,oi ' 2 


2j v n 


(A.16) 


Lj (ityx y 


</^U 4 

y/kn 


<v^ 

VK 


^n 2( t>n + ^^iaxE ra . Vn {U 2jVn+ (2 l -l)A n ~ E'K))' 


< a: 


Similarly, by the property of Ut(u n ), (6.22) and (6.29) in Jacod and Todorov 
(2014), and Holder’s inequality, we have 


„ Ef=r'co0',o / r (i/(fc n -i))Eitr 1 (ci(j,o-CiO‘,o) 

—-I-—— - 


2jvn 


(^n 1) Uj (Vjn) V 


' fan 1 ^ 
Uj (u r j 


< 


(A.17) 


Efer^^cgao 

M (/c n - 1) 2 E 2 K) 


x W Ejr, 


(i/(^-i))Ef=r 1 ic^(j,o-Co(j,or 2 




Uj (Un ) 


< A- 




^/kyi 

Equations (A.16) and (A.17) yield 

■^ 2 . jv n ^0,j{ u n)^,l,j(u n ) 

(A.18) 

„ (V(tn - i^Efer’coaoti/fe - iDEferUiQ-.o 

^ 2i ’" Uj(u n )Uj(u n ) 


+ r 


Ji 
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where rj satisfies \rj\ < Ky/u^u\ / \/^n■ By the dehnition of we have 


(i/(t» -1)) Efcr 1 Coo,o(i/(<v, -1)) ESr 1 cib'.o 

Uj (lL n ) Uj (lL n ) 

kn 1 


k n 1 


E ?2r 


2jvn 


^ n 

(A ' 19) = (fc„-l) 2 ^(«n)^M U 


kn—2 


+ X/ ^2 J -«„Co(j»OCiO'^ + 1 )- 


1=1 


By Lemmas 11-12 in Jacod and Todorov (2014), we have 

^■^CoO’.OCiO'.O 


= E^,„Co(i, i) 


(A.20) 


= E T 2jVn cos 


X cos u. 


(u n ^ 


+2Z+1 J 


v n . 

2jkn+2l^ 


A 1/2 

\ AA n 

An y_ An y 

AA 2j/c n +2/ vV ^2j/c n +2/-l^ 

! A 172 


Uj ( U n ) Uj ( Un ) ~\~ r2 j j 

where r 2 j satisfies \r 2 j\ < Ku^y/v^- Since cos(x)cos(y) = ^(cos(x + y) + 
cos(x — y)), we have by Lemma A.4, 


E Tly„ n cos ( U ' 


\n y— A n y 

^2.jfc n +2Z+l ^ 2jk n +2l^ 


A. 


1/2 


/ A n y _ A n y 

/ ^ A 2ifc„+2Z y ^ ^2?fc„+2Z-l^ 
COS «„— - 773 -^- 

V Ai /2 


= -£V„. cos 

O 2j-j n I 


/ A n y A n y 

/ ^2jfc n +2Z+l^ ^ A 2ifc„+2Z-I ys - 

1 Un 1/2 “ ?t a 1/2 




A 


(A.21) 


An y A n y A n y 

, ^2jk n +2l+l^ ^2jfc„+2i A ^2jk n +2l-l^ 

I LUO | U n . £ LLfi , /o L" 0 - 77 , 


A 


1/2 


A. 


1/2 


A 


1/2 


= o E T 2jVn ( ex P (~ulc 2 jv n +( 2 l- 2 )A n ) + exp {-3u 2 n c 2jVn+( 2i-2)A n )) + r 3J 


= 7 (exp (~u 2 n c 2jv J + exp (-3 u 2 n c 2jVn )) + r 3J + r 4 j, 

where |r, 3 j| < KAh and |r 4 j| < Kv? n v n by second-order Taylor expan¬ 
sion on e x for x around the origin and (S.1.3) with V = c. Now substituting 
(A.21) back into (A.20), we have 

E r 2j vn Co (A 0Ci(j, 0 = | (exp (- u 2 n c 2jVn ) + exp (-3 u 2 n c 2jVn )) 
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(A.22) 

- exp {— 2 u n C 2 jVn ) + , 

where |r5^ | < K(y/v^+ A, 1 , /3 ^ 2 ). Similarly, we have 

E T 2jVn Co(j,l)Ci{jJ + 1) 

(A.23) = i(exp {-u 2 n c 2jVn ) + exp (-3 u 2 n c 2jVn )) 

— exp (— 2 u^ l C 2 j Vn ) + tqj, 

where \rej\< K{y/v^ + An /3 ^ 2 ). Substitute (A.22) and (A.23) into (A.19), 
and then substitute the latter into (A. 18), and we have 


E J r 2 :j Vn ^0,j{ u n)^l,j( u n) 


(A.24) 


exp (-ulc 2jVn ) + exp (-3 u 2 n c 2 j Vri ) - 2exp (-2ulc 2jVn ) 
(k n - 1 )Uf(u n ) 

u n4 jVn + r 7j , _* _ < c ljv n + r S j , 

1 'n 7 1 'j ? 


1 * 

+ r,- 


(/c n - l)C/|(u n ) J fcn -1 

where |r 7j -1V |r 8 j | < KvP n , |r*| < |r,- | + 1 r2j 1 + 1 r3j ' 1 1 r5j 1+1 r&j 1 + ■ Now 

a combination of (A.24) and (A. 14) yields 

4, v 2 lt/( 2 v„)\-l 

E E J~ 2 jv n (Co ,j (u n ) - E J~2jv n Co ,j (u n )) 


A n n n 


(A.25) 


i=o 


X (Cl ,j( u n) E J r 2 jv n ^,j( U n)) ~ 2 f C s ds. 

■JO 


This proves the third equation in (A. 12). 

By Lemma A.2, we also have 

[t/( 2 Vn )]-l 

E \xf\>e) 

l=o 


1 [t/(2«„)]-l 
— 7 ^^'DnlX 


n,fc|3 


1=0 

~ [t/(2«„)]-l 
l=o 


2 Vn 


A i /2 y «n 


-.6 Ejr 2jvn ICfej'^n)! —t 0. 


This proves the Linderberg condition [equation four in (A. 12)]. 
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Taking k = 2 and = 1 in Lemma 15 of Jacod and Todorov (2014), we 
have 

[t/(2u n )]-l 

Ejr *jvu Xj' k ( M 2(j+1)** - M 2 jVn ) 

3= 0 

\t/(2v n )\—\ ^ 

= zL ~ 1/2 & J( U J ) ( M 2 (j+1)«„ - 0. 

j=0 ^riAn 

This proves the final equation in (A. 12) and completes the proof of the 
bivariate central limit theorem with stable convergence. □ 


Proof of Theorem 3.2. Let T nX = C 0 (u n ) - A 0 (u n )^ -C t - (Ci(u n ) - 
Mun)?~C t ). By (2.5), 


(A.26) 


T n = 


T n i - TnA^ 2 


ini 


c t + MOt+o P {ti! 2 ) c t + o P (i) 


+°p(Ay 2 ). 


Then Theorem 3.2 is a straightforward consequence of Theorem 3.1, (A.26), 
the stable convergence mode and the continuous mapping theorem. □ 


Proof of Theorem 3.3. By Theorem 3.1, C\{u n ) = C t + A 0 (u n )™ + 

Op(A n ' ) = Ct + o p ( 1). This shows that the denominator of kt converges to 
in probability. By (A.2), we have 


(A.27) 
where 


c k Au n ) - 


(sinh (ultf(u n ))f 


3X '~ n ' {k n -l)ul 

= c 2jv n + ^j,l( u n) + Cj, 2 ( u n ) + C^ 3 (u n ), 


k _ (£k,j( u n) 


C 3 ,1 = 


C 7 -2 — 2C2j 




+ 


£kj(“n) (sinh(u2 c ^(u n ))) : 


2u2 


(fc n - 1 )u 2 n 


+ (^fcj(Wn)) 2 , 


Cfcj(un) 4jK) (sinh(<cj(u n )))^ 

O I 


J ’ 2 1 2u£ (A; n -1K 

Cj, 3 ( u n) =4c2 jVn uP- 2 A 1 n -P/ 2 a 2jVn + (2^ _2 A 1 -P /2 a 2 jVn ) 2 
By (A.l), 




[t/(2v n )]-l [t/(2«n)]-l 

. 7=0 j =0 
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By Lemma A.2, o 


[t/( 2v n )\-l, 


2v n {^±) 2 = o p { 1). On Q{k,n,t), - 

“'ll " “n 


( U «) 


(sinh(« 2 cH« n ))) 2 , . , ,, zs I 2 i, 

— - 1 is bounded by A/u n , hence 


(fc„-l)u 2 


Cfej(^n) (sinh(u2 c ^(u n ))) 
2 ul 


(kn ~ 1 )«n 




K 

Zk,j( u n) 

(sinh {u 2 n c k {u n ))) 2 

7/2 

LL n 


{k n - 1 )u 2 


/A x K 

(A.28) <-j 




LQ(k,n,t) 


*k,j 


,U r 


£k,j [ 


U r 


2ul Ej ^ 2ul 


+ 


E ?2t 


tljiO U j (2u n ) + l-2Uf(ul) 


2jv n 


2u 2 


4u 2 (k n - 1 )Uf(u r 


+ 


Uj(2u n ) + 1 - 2 U 2 {u 2 n ) (sinh {u 2 n c k j{u n ))Y 


4«2 (L n _ l)u 2 {u r , 


{kn - l)u 2 


*r2(/c,n,t) • 


By the property of Uj{u n ) and the definition of c k {u n ), the expectation of the 
third absolute value conditional on jv n is smaller than K {un A n ' /k n + 
u~ 2 /kfl 2 + A}/ 2 (f> n /k n ). By Lemma A.2, the second absolute value is smaller 
than Ku^An ' 2 (f) n . By Holder’s inequality and Lemma A.2 with q = 4, the 
expectation of the first absolute value conditional on J-2jv n is smaller than 
K{un/k n + y/v^). In summary, we conclude that 


[t/(2u n )]-l 

(A.29) X] 

3=0 



(sinh {u 2 n c k {u n ))) 2 
k n u 2 


2 


° p ( !)• 


By (A- 1 ), Y}j=Q Vn)] 1 ‘ lv n{r k ,j{u n )) 2 I U c( <Knj f ) — Op (1)• On Q{k,n,t), \r kjj \ < 
K \ikjun)\ ' gy Lemma A.2 with q = 6, we have Xo=o 2 ° n ^ _1 2v n {r k j{u n )) 2 x 
I(l(k,n,t) = o p (l). Combining all the results of the terms on the right-hand side 
of the decomposition of c k x , we have Y!j=Q Vn ^ - * 2v n c k , = o p (l). Similarly, 
one easily proves that Y^jL 1 2v n c k 2 = o p (l). By boundedness of c and 
a, Y^j=Q Vn ^ 1 ^nCj i 3 = o p (l). This shows that 


[T/(2v n )\ 1 2 ^ 

(A.30) i nk = V 2v n cl jVn + o p {l) = c 2 ds + o p { 1), fc = 0,l. 

j=o • /o 

This shows that the numerator of Ky converges to 4 J^ c 2 ds in probability, 
and hence kt itself converges to kt in probability. On the other hand, by 
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1 /2 

Theorem 3.2, T n /A n converges to Gt stably. By the stable convergence 
mode, T n converges to standard normal distribution stably. □ 


A.2. Proof of results under H\. In the sequel we assume that A is a 
pure-jump process. We rewrite 

log Uj(u n ) log(l + 6j(«n)) 


Cj(u n ) = -- 


(A.31) 


3 V 2 

tlrr, 


Ut 


-2u /3 - 2 A 1 - /3/2 a9- - +f, ■ 

— z,u n. L - X rt. u 2jVn n 2 ' k,]i 


ui 


k = 0,1, 


where < K^ j{u n ) /n 2 on Q(k,n,t). Recall the definition of 
(A.26), and we have 


m 


(A.32) 
where 


[t/(2u»)]-l 

T n i = 2n n ^ 

3=0 




+ A; 


Ut 


■n,t j 


[i/(2u„)]-l 


Rn,t — 2n n ^ ] 

3=0 

(A.33) 


(R)j - rij) 


+ 


(sinh(u 2 c)(u n ))) 2 - (sinh(u 2 cRu n ))) 


U 2 n (k n - 1) 

Similar to Lemma A.2, we have the following. The proof is provided in 
the supplementary material [Kong, Liu and Jing (2015)]. 


Lemma A. 5. Assume Assumptions 3.1, 3.3 and A.l, and suppose u n is 
bounded, so we have on the set {Ct = 0}, 

I ^T 2 j Vn €k,j(u n ) | 

(A.34) < r / 2 ) A (( 3_/3 )/ 2_£, ) A ( 1_ ^/( 2 ( /3+1 — r ))) 

-j/2_ ^ il a 1/2 

and if further k n A n —>• oo /or any e > 0, and sup n n ~t n < oo is satisfied, 


(A.35) 


E ?Vv n &j( u n) < K 




The following lemma gives the convergence rate of the terms on the right- 
hand side of (A.32). The proof can be found in the supplementary material 
[Kong, Liu and Jing (2015)] to this paper. 
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Lemma A.6. Assume Assumption 3.1, 3.3 and A.l, and suppose u n is 
1 /2 

bounded and k n An —)• 0, so we have on the set {Ct = 0}, 


(1) 


(A.36) 


[t/{2v „)\-1 

2 Vn Ejr 2 jVn 

3=0 


Co j'(^n) Cl,j(^n) 


Ut 


< j R' u -2( A (l-r/2)A((3-^)/2- £ ')A(l-^/(2(/3+l-r))) + u ^3/2-/3/2). 


( 2 ) 


(A.37) 


[t/(2v„)]-l 

^ ! ^^2jv„ 

3=0 


^0 ,j( u n) Cl ,j{ u n) ^J r 2 jv n ^ 0 ,j{ u n) Cl ,j( u n)) 

ul/{2v n ) 


< Ku-\uiA 2 -^ 2 + A Cf—/ 2 ) A C(5-/3)/2-eO ); 


/or any e' > 0; 

(3) 

(A.38) 


ddnt — Op 


ui^A^ 2 


Proof of Theorem 3.4. We first prove the first equation. By (A.31), 
we have 


(A.39) 


Co (n n ) C i(u n ) — Tn,i 


[t/(2v n )] — l 

2 V Co ,j\ u n) Si ,j\ u n) 

3=0 


Ut 


+ Rn,t- 


Now by Lemma A.6, we have 

(A.40) Co(u n ) — C\(u n ) = Op(5n t 1 + yj 5n t 2 + bn, 3 ), 

where 

S n x = u -2^ A a-r/2)A((3-/3)/2- £ ')A(l-/3/(2(/3+l-r))) + u P/±V2-/3/2^ 

ysu= 

4 .,3 = <‘L 2 Al _|3/2 /t„. 

Now, notice that: (1) 1 — r/4 > 1 — r/2 > 1 — /3/2(/3 + 1 — r); (2) 5 ni 3 > 
n£- 2 A^ /2 ; (3) ^ (4) ^ /2 _ f , < v&/{k n ti! 2 ~ e ') < *4 

U n A n ' 
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-2 A 3 / 1 - 0 / 2 -e' A—R /9 i A 1 / 2 , 

(5) M ^/t- 2 A i-/ 3/4 = Un n ui - J/ 4 + e ' < -K'- By choosing e' > 0 small 


n k n A„ 

enough and the conditions on u n and k n , we have 


£•„(»„) - ft ft) = o„(»- 2 a;,-«w«--» + u^A^-KR). 

Next, we prove the second equation. By (A.31), we have 


(A.41) 


(< A («„)) 2 - ( 2 u%~ 2 Al l -P^a 2 j Vn ) 2 - ( M ? A ) 2 - ft ,,) 2 

v \ u n / 


= 2(2ut 2 ^ n - 0/2 a 2jVn )(- 


£k, j(Un) ~ 
,2 “C 




-2 


£fc,7 ( w n) 




r fc, 7 ‘ 


Now we use several steps to show that under Hi the principal term of 
(i c k Au n )) 2 is (2u^ 2 Al l ~ l3/2 a2jv„) 2 and Ejr (tf(u n )) 2 < Kv^~ A Al~^. By 


Lemma A.5, we have 

(A.42) (LiLL) 2 < Kiuf 2 AT» 2 ) 2 


kn A^/ 2 Un 13 


u: 


n fc 2 A n 


3/2—/3/2 : 


which is o p ((un 2 An ^/ 2 ) 2 ) by the conditions on k n and u n given in The¬ 
orem 3.4. By Lemma A.5 and (A.42), we have on £l(k,n,t) (on which 
I'ffcjl < K^lj/ul and is bounded), 

(rfcj) Ifi(k,n,t) — KEjr 2 jv n \fk,j\lQ(k,n,t) 

£k,j{u n ) 


(A.43) 


< 


2jvn 




< K(ut 2 Ai-KR) 


■>kaT 4 0 

< KA'J 2 - 0 ' 2 


By (A.42) and (A.43), we have by Holder’s inequality, 




2jv n 


(A.44) 


{ 2ut 2 Ai-^a 2 , Vn ) - 


£,k,j {Un) 


< K(u 0 - 2 A 2 -^r 


Ut 


4-/S/2 

U n 


+ r k . 


iQ(k,n,t) 


+ ■ 


Ut, 


,3/2 A l-/3/4 , 2 a 1/2 I 


and 

Ej~2jvn 
(A.45) 


Ck,j(u n ) 


Ut, 


r k ,j | Iil(k,n,t,) — ■EEj r 2 j v 


Ck,j {Ur, 


< K«~ 2 A 2 - 0 ' 2 ) 


2 k n A n 


1/2 


4-/3 

U n 


< k 2 aI ,2 - 0/2 ' 
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Combining (A.42)-(A.45) yields that 

(A.46) E T21 J(<*(u n )f - {2<- 2 A;-'>/ 2 ) 2 |= 0(1), 

where o(l) holds uniformly in j. 

By the form of Cj(u n ), we have u 2 |Cj (u n )\I^ k n ^ < K, and hence by 
Taylor expansion on the exponential function, we have 

(A.47) (sinh (ultf (u n ))) 2 In( k ,n,t) < W)%(fc,n,t) < K. 

By virtue of (A.47), we have 

'(sinh(^cj'(^n))) 2 /o(fc,n,t) N 2 


E ^2r 


2jvn 


u n{kn ~ 1 ) 


(A.48) 


< 


u : 


K 

4 ; 2 ^^2 jv n (sinh(?i n Cj 
n K n 


< §E^(^(„ n )) 2 I„f iin:t) < A(„6- 2 AyO/ 2 ) 2 , 


1,2 

r^n 


and further by the Cauchy inequality, 
(sinh {u 2 n c k Au n ))) 2 


(A.49) E T , 


cAu n )- 


K 


*>(*,„,*) < ^{vt 2 K~ P/2 ) 2 - 


3 ul(k n - 1) 

Now combining (A.46), (A.48), (A.49) and (A.l), we have 

/[</( 2 v n )\-l \ 

in,k = (2uf” 2 A^- /3/2 ) 2 ^ Y2 2v nal jVn + Op( 1) j 

= (2^ -2 A^ /3/2 ) 2 ^ a 2 ds + o p (l)^ , 
for k = 0,1. This proves the second equation of Theorem 3.4. □ 


Proof of Corollary 3.1. Part 1 is a straight consequence of Theo¬ 
rem 3.3. Now we prove part 2. By Theorem 3.4, we have by the condition 
on 7 „, 


-7n + o r K 2 + uiP-Gy-w 4 ) 

4 u%~ 2 a 1 t 0/2 \Jf*a% ds + Op( 1) 

7n(l + °p(l)) 

4u^“ 2 A^, _/3/2 y // /qO 2 ds + Op(l) 


(A.50) 
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Since 7 n u 2 ^ 2 A ^ 4 1 ^ 2 —> 00 and 


2-/3/2 a /3/4-1/2 
Jj n __ 

2-/3 a /3/2-1 


< » 2 / 2 a 7-« 4 -» 0 , 


7n(l + °p(l)) 


4un 2 An ^ /2 ^/ Jq a 2 s ds + Op(l) 


— 00 . 


This proves part 2 on the performance of the power of the test. □ 


SUPPLEMENTARY MATERIAL 

Supplement to “Testing for pure-jump processes for high-frequency data” 

(DOI: 10.1214/14-AOS1298SUPP; .pdf). This supplement contains technical 
proofs of the Lemmas A.4-A.6 as well as some interesting supplemental 
lemmas. 
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